Metal-insulator transitions in cyclotron resonance of periodic nanostructures due to 

avoided band crossings 
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A recently found metal-insulator transition in a model for cyclotron resonance in a two-dimensional 
periodic potential is investigated by means of spectral properties of the time evolution operator. 
The previously found dynamical signatures of the transition are explained in terms of avoided band 
crossings due to the change of the external electric field. The occurrence of a cross-like transport is 
predicted and numerically confirmed. 
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I. INTRODUCTION 

Since the pioneering work of BlochEl and LandarJl the 
investigation of spectral properties and. transport fea- 
tures of Bloch electrons in magnetic fieldaJ have attracted 
much attention. While infinitely degenerate Landau lev- 
els are found for free electrons in a static magnetic field, 
a spatially periodic potential without field yields a band 
structure in the spectrum. The general situation of a pe- 
riodic potential and-a, magnetic field is still a matter of 
intensive researchoT]. A particular simple realization for 
a two-dimensional (2D) periodic potential is, e.g., given 

by 



V{r) = Vx cos ( ^) + Vy cos 



(1) 



where a, b are the periods in the corresponding di- 
rections. In a one-band approximation this potential 
leads to the well-known Harper modelQ, where a metal- 
insulator transition (MIT) is observed upon changing the 
modulation amplitudes from Vy > Vx to Vx > VyU. At 
the critical point, i.e., Vx = Vn the spectrum and the 
eigenfunctions- are , rn ultifractalsll^ and give rise to anoma- 
lous diffusioiJB'BEJ. The one-band Harper model, how- 
ever, has an integrable classical limit, whereas Bloch elec- 
trons in a magnetic field show chaotic dynamics in the 
classical limit. In quantum mechanics this is refiected in 
the coupling of Landau bands. This coupling leads to 
avoided band crossings, which induce an unusual trans- 
port phenomenon: For Vx ^ Vy electrons show ballistic 
transport in the direction of strong potential modulation 
and localization in the direction of weaker modulationO. 

Analogous phenomena can be found in driven systems 
with a chaotic classical limit even within the one-band 
approximation. A welL studied example is the kicked 
Harper model (KHM)L3ilj. For small kicking -strength 
its behavior is analogous to the Harper modelEj, while 
for increasing kicking strength avoided band crossings 
due tp-.the classical non-integrability induce a variety of 
MITsEfl. Recently, it was demonstrated that a variety of 
kicked models, including the KHM, can be realized for 
electrons on a lattice in the presence of a magnetic field 



driven by a smooth electric field@. The effective kicking 
is a result of a resonant interaction between the electronic 
motion and the driving field. These models correspond 
to cyclotron reso nanc e experiments in antidot arrays and 
in organic metalgc30. By investigating the dynamics, a 
localization-delocalization transition induced by the am- 
plitude of the electric field was obscrvcdEil and remained 
unexplained so far. 

In this paper we investigate these dynamical proper- 
ties in terms of spectral properties in a wide range of 
driving strengths. We show that the previously found 
localization-delocalization transition is a consequence of 
avoided band crossings due to the change of the driving 
field strength. Moreover, the spectral analysis predicts 
a cross-like ballistic transport, which is numerically con- 
firmed. We start in Section II by briefly introducing the 
model and explain the numerical approach in Section HI. 
In Section IV we present our main results, which give a 
clear picture for the MIT in terms of the spectrum. We 
finally summarize our findings in Section V. 



II. THE MODEL 

The one particle Hamiltonian, describing the cyclotron 
resonance of an electron in the 2D periodic potential V(v) 
of Eq. (|^) can be written in the following form: 



^= :7^(p-eA(r,i))' + V(r) , 
2m* 



(2) 



where p = {px,Py) is the momentum of an electron 
with effective mass m* and charge e moving in the xy- 
plane perpendicular to the constant magnetic field B 
pointing in z-direction. We choose the Landau gauge, 
which yields a vector potential A(r, t) = (^ sin i/t, xB — 

cos vt) containing also the alternating electric field 
E — {~Ex cos i/t, Ey sin i^t) with frequency v. 

Applying a gauge transformation to the Hamilto- 
nian (^ , the time dependence-can be shifted to the peri- 
odic potential V{r) V{r,tjB: If the fields B, Ex, Ey 
satisfy uj*Ey = vEx then one can rewrite the Hamilto- 
nian as H = Ho(p,a;) -I- V{r,t), where 1-Lq = Ho{p,x) 
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describes the cyclotron motion with frequency u* = 
eB/m*. The perturbation V{r,t) splits the infinitely de- 
generate Landau levels of Ho into a mini-band structure. 
For strong electric fields ap4rif = f this Hamiltonian H 
reduces to a kicked systemc3€J 

H = Lcosp + K { cos((7 - Ko)S{t - 2n + 1) 

71— — OQ 

+ cos(g + Ko)d{t - 2n) } , (3) 

where the canonical operators q and p corresponding to 
the spatial coordinates x, y satisfy [p, q] — ~i2nhcs and 
kq = {2neEy/ {hm*) — 7r/4) mod 2tt\ i, K denote the 
rescaled amplitudes of the potential (|l|). The effective 
Planck constant Tioff = Ti/(aheB), corresponds to the in- 
verse number of magnetic flux quanta through a unit 
cell a X 6. The_aystem described by Eq. (^) has a chaotic 
classical limited and for kq = reduces to the KHM. 



III. NUMERICAL APPROACH 

The time evolution operator for one period correspond- 
ing to (H) reads 

;y ( \ ( I ^\ 

LA ~ exp — - — cosp exp — - — cos(g — kq) 

V "off / V 'ioff / 

( \ ( iK . A 

X exp — - — cosp exp — - — cosyq + kq) . (4) 

V "off / V "cff / 

In order to obtain its spectrum, one has to approximate 
the irrational effective Planck constant by a rational ap- 
proximant, ?iefT = 2TrMJM. Following standard tech- 
niques for kicked system£3, this leads to a x A^-matrix 
for the evolution operator U, which depends on two Bloch 
phases Oq,0p £ [0, 2tt]. These Bloch phases arise from pe- 
riodic boundary conditions in the q- and p-direction, re- 
spectively. Diagonalization oiU for all pairs {9q, Op) gives 
its spectrum. Keeping 9q fixed and varying Op yields a 
part of this spectrum, from now on referred to as the 
"p-spectrum" . Eigenfunctions, which are localized in p- 
direction will be associated with very narrow bands in the 
p-spectrum, whereas extended states give rise to broad 
bands. Since for higher approximants of Tics the width 
of the narrow bands decreases exponentially, we will call 
them levels. In an analogous way one can define the "(7- 
spectrum" , which reflects the properties of the eigenfunc- 
tions in g-direction. Figure ^ shows a q- and p-spectrum, 
where one observes broad bands and levels. 

In order to investigate dynamics, we will study the 
spreading of wave packets, initially (5-localized in either 
q- or p-direction by determining the corresponding vari- 
ances Mqit) and Mp{t). For the spreading in p-direction 
we use periodic boundary conditions in ^-direction after 
one unit cell, giving a discrete lattice in momentum space 
with spacing ?iefT- An analogous construction is used for 
the spreading in g-direction. When considering transport 



in phase space in both directions, we use periodic bound- 
ary conditions after r unit cells in g-direction, which gives 
r different Bloch phases 9q — 2Tr j/r (j = 0, . . . , r — I). 
The complete wave packet is then obtained by combining 
the wave packets for each Bloch phase and it is defined 
on a grid with spacing hcs/r in momentum space. 

IV. RESULTS 

We will study the spectrum of Eq. as a function 
of Kq and its consequences for dynamical properties. b| 
particular, we will focus on the parameters used in Ref.tj 
and show, that they reflect the generic behavior of the 
model (|L, 

In RefEil for X = L = bandhcs ^ 2tt/{7 + ag), where 
(Tg = (\/5 — l)/2 denotes the golden mean, a transition 
from diffusive to ballistic dynamics was foimd, when kq 
was varied from to iran (cf. Fig. 1 of Ref.El). In Fig. |l| 
we present the q- and p-spectrum for one period of kq. 
Apparently, the duality of q- and p-spectrum is conserved 
for all Kq, i.e., the g-spectrum shows wide bands, where 
the p-spectrum is level-like and vice versa. 

For Kg = the variances increase almost linearly 
in both q- and p-direction, corresponding to diffusive 
spreading in either direction (top of Fig. ^). This 
well-known behavior of the KHM is a direct conse- 
quence of the- miiltifractal nature of spectrum and 
eigenfunct ionsE3i]£l. 

For Kq = TTCTg bauds occur in the p-spectrum explaining 
the ballistic transport in p-direction observed in Ref.tiJ. 
There are, however, also levels, which are related to 
bands in the g-spectrum, predicting a ballistic spread- 
ing in g-direction. This is confirmed in Fig. |^. All 
eigenstates are extended in either of the directions and 
localized in the other one. A typical initial wave packet 
will excite both types of eigenfunctions. The part of the 
wave packet consiting of eigenfunctions extended in the q- 
direction gives ballistic transport in the g-direction. The 
other part, which is a superposition of eigenfunctions ex- 
tended in the p-direction, yields ballistic transport in the 
p-direction. Therefore one finds a superposition of bal- 
listic transport in q- and p-direction, namely a cross-like 
transport (bottom of Fig. |pD). The cross- like form of the 
wave packet is clearly different from the isotropic form 
for Kq = (Fig. ||a). There are several other trans- 
port features one can observe from the Husimi plot in 
Fig. The weight of the wave packet spreading bal- 
listically along the p-direction is much bigger than the 
corresponding weight for the g-direction. This is a signa- 
ture of the spectrum, namely that there are more bands 
in the p-spectrum than there are bands in the g-spectrum 
at Kq — TTCTg, as can be seen in Fig. ^ It should be noted, 
that the velocities in p- and q-direction (which can be in- 
ferred from the variances Mp(t) and Mq{t)) are in general 
not related to the number of bands but are determined 
by the band widths. 
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Now we will turn away from the symmetry line K = L 
and focus on the parameters ivT = 2, L = 4, ?icff — 
27r/(7 + (Tg). In Fig. |^ we show the q- and p-spectrum 
as a function of kqj which again appears to be dual. For 
the KHM (ko = 0) the g-spectrum shows broad bands, 
whereas the p-spectrum is level-like, indicating-jbransport 
in g-direction, but localization in p-directionlU. This is 
confirmed by the dynamics (Fig. IJa). Variation of kq in- 
duces avoided band crossings leading to the occurrence 
of bands in the p-spectrum (see magnification in Fig. 3). 
This is associated with a localization-delpcalization tran- 
sition in p-direction as observed in Ref.cll. As there are 
still bands in the g-spectrum one again finds the generic 
cross-like transport (Fig. ||b). Here, the weight of the 
wave packet spreading along the g-dircction is bigger than 
the corresponding weight in p-direction, as is implied 
by the much larger number of bands in the q-spectrum 
(Fig. H) . Another characteristic of the wave packet is the 
increased localization length in p-direction of the part 
of the wave packet spreading in g-direction compared to 
Fig. This is due to the increased locaHzation length 
of the eigenstates corresponding to the levels in the p- 
spectrum. This change in the localixajbion length can be 
explained by avoided band crossingaia. 

V. SUMMARY 

In conclusion, we have shown here, that the previously 
found MIT in cyclotron resonance in a 2D-periodic poten- 
tial in the presence of a magnetic field can be fully under- 
stood by investigating the spectrum of the corresponding 
time evolution operator. We have presented numerical 
evidence, that the duality of the spectra is conserved un- 
der changes of the driving strength. We|-i*fere able to 
explain the previously observed dynamica22l, which has 
shown a transition from diffusive to ballistic transport 
as well as a transition from localization to ballistic trans- 
port in the p-direction. Furthermore, the spectra allowed 
us to predict the transport behavior in g-direction lead- 
ing to a cross- like transport. These predictions were also 
confirmed numerically. We would like to mention, that 
cross-like transport is generic for systems with dual spec- 
tra and a chaotic classical limit. 

The observed transport phenomena in phase space of 
the effective Hamiltonian (^) correspond to analogous 
transport phenomena in the xy-plane. The required con- 
ditions for realization of the observed MIT m^ well be 
achieved in experiments on cyclotron resonanceu of a 2D 
electron gas embedded in iateral superlattices fabricated 
on GaAs heterostructuresEJ. 
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FIG. 1. The p-spectrum (top) and g-spectrum (bottom) of 
the evolution operator (^ for K — L — 5 vs. kq (one period 
[—77/2,71/2] is shown). The approximant for the irrational 
effective Planck constant is fteff = 27r- 34/259. Arrows denote 
the values kq = and kq = TrcTg used in Fig. ^. For kq = 
Eq. (^) reduces to the KHM and the spectrum is multifractal, 
while for kq 7^ it consists of bands and levels. The spectra 
appear to be dual to each other, since broad bands in one of 
them are associated with levels in the other. 
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FIG. 2. (a) Variances Mq{t), Mp{t) and the Husimi plot 
of a 2D wave packet after 100 kicks (one point per unit 
cell and r — 89 cells in each direction) for K — L = 5, 
fioff ~ 27r/(7 + (Tg), Ko = 0. (b) same for kq — nag. In the 
first case, the variances show a diffusive-like behavior while 
the wave packet spreads ballistically in both directions in the 
second case. This is also reflected in the Husimi plots, where 
one observes that the wave packet spreads isotropically in 
phase space in (a) and shows a cross-like spreading in (b) . 



FIG. 3. Same as in Fig. |l] for if = 2, L = 4. For kq = 
the p-spectrum consists of levels only and the g-spectrum 
shows bands. Magnifications [hctt = 27r55/419) show part 
of the spectra around ko ~ nag, where one finds bands in the 
p-spectrum and levels in the g-spectrum. Duality again seems 
to be conserved. 



FIG. 4. Same as in Fig. | for if = 2, L = 4. (a) Ballistic 
spreading in g-direction and localization in p-direction is ob- 
served, (b) Ballistic spreading in both directions leads to a 
cross-like transport. 
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